The emerging complexity of Open Science: assessing Intelligent Data Openness in Genomic Anthropology and Human Genomics.

In recent decades, the scientific community has become aware of the importance of science being effectively open in order to speed up scientific and technological progress. In this context, the achievement of a robust, effective and responsible form of data sharing is now widely acknowledged as a fundamental part of the research process. The production and resolution of human genomic data has steadily increased in recent years, mainly due to technological advances and decreasing costs of DNA genotyping and sequencing. There is, however, a downside to this process due to the huge increase in the complexity of the data and related metadata. This means it is advisable to go beyond traditional forms of sharing analysis, which have focused on data availability only. Here we present a pilot study that aims to complement a survey on the availability of data related to peer-reviewed publications with an analysis of their findability, accessibility, useability and assessability (according to the "intelligent data openness" scheme). Sharing rates in genomic anthropology (73.0%) were found to be higher than human genomics (32.4%), but lower than closely related research fields (from 96.8% to 79.2% for paleogenetics and evolutionary genetics, respectively). We discuss the privacy and methodological issues that could be linked to this finding. Comparisons of sharing rates across a wide range of disciplines has suggested that the idea of human genomics as a forerunner for the open data movement should be questioned. Finally, both in genomic anthropology and human genomics, findability and useability were found to be compliant with the expectations of an intelligent data openness, whereas only a minor part of studies met the need to make the data completely assessable.

Berbers and Arabs: Tracing the genetic diversity and history of Southern Tunisia through genome wide analysis.

OBJECTIVES: Tunisia has been a crossroads for people from Africa, Europe, and the Middle East since prehistoric times. At present, it is inhabited by two main ethnic groups, Arabs and Berbers, and several minorities. This study aims to advance knowledge regarding their genetic structure using new population samplings and a genome-wide approach. MATERIALS AND METHODS: We investigated genomic variation, estimated ancestry components and dated admixture events in three Berber and two Arab populations from Southern Tunisia, mining a dataset including Middle Eastern, sub-Saharan, and European populations. RESULTS: Differences in the proportion of North African, Arabian, and European ancestries and the varying impact of admixture and isolation determined significant heterogeneity in the genetic structure of Southern Tunisian populations. Admixture time estimates show a multilayer pattern of admixture events, involving both ethno-linguistic groups, which started around the mid XI century and lasted for nearly five centuries. DISCUSSION: Our study provides evidence that the relationships between genetic and cultural diversity of old and new inhabitants of North Africa in southern Tunisia follow different patterns. The Berbers seem to have preserved a significant part of their common genomic heritage despite Islamization, Arab cultural influence, and linguistic diversity. Compared to Morocco and Algeria, southern Tunisian Arabs have retained a higher level of Arabian ancestry. This is more evident in the semi-nomad R'Baya, who have kept their original Bedouin lifestyle, than in the population from Douz, who have undergone multiple events of stratification and admixture.

Inter-individual genomic heterogeneity within European population isolates.

A number of studies carried out since the early '70s has investigated the effects of isolation on genetic variation within and among human populations in diverse geographical contexts. However, no extensive analysis has been carried out on the heterogeneity among genomes within isolated populations. This issue is worth exploring since events of recent admixture and/or subdivision could potentially disrupt the genetic homogeneity which is to be expected when isolation is prolonged and constant over time. Here, we analyze literature data relative to 87,815 autosomal single-nucleotide polymorphisms, which were obtained from a total of 28 European populations. Our results challenge the traditional paradigm of population isolates as structured as genetically (and genomically) uniform entities. In fact, focusing on the distribution of variance of intra-population diversity measures across individuals, we show that the inter-individual heterogeneity of isolated populations is at least comparable to the open ones. More in particular, three small and highly inbred isolates (Sappada, Sauris and Timau in Northeastern Italy) were found to be characterized by levels of inter-individual heterogeneity largely exceeding that of all other populations, possibly due to relatively recent events of genetic introgression. Finally, we propose a way to monitor the effects of inter-individual heterogeneity in disease-gene association studies.

Overcoming the dichotomy between open and isolated populations using genomic data from a large European dataset.

Human populations are often dichotomized into "isolated" and "open" categories using cultural and/or geographical barriers to gene flow as differential criteria. Although widespread, the use of these alternative categories could obscure further heterogeneity due to inter-population differences in effective size, growth rate, and timing or amount of gene flow. We compared intra and inter-population variation measures combining novel and literature data relative to 87,818 autosomal SNPs in 14 open populations and 10 geographic and/or linguistic European isolates. Patterns of intra-population diversity were found to vary considerably more among isolates, probably due to differential levels of drift and inbreeding. The relatively large effective size estimated for some population isolates challenges the generalized view that they originate from small founding groups. Principal component scores based on measures of intra-population variation of isolated and open populations were found to be distributed along a continuum, with an area of intersection between the two groups. Patterns of inter-population diversity were even closer, as we were able to detect some differences between population groups only for a few multidimensional scaling dimensions. Therefore, different lines of evidence suggest that dichotomizing human populations into open and isolated groups fails to capture the actual relations among their genomic features.

When data sharing gets close to 100%: what human paleogenetics can teach the open science movement.

This study analyzes data sharing regarding mitochondrial, Y chromosomal and autosomal polymorphisms in a total of 162 papers on ancient human DNA published between 1988 and 2013. The estimated sharing rate was not far from totality (97.6% ± 2.1%) and substantially higher than observed in other fields of genetic research (evolutionary, medical and forensic genetics). Both a questionnaire-based survey and the examination of Journals' editorial policies suggest that this high sharing rate cannot be simply explained by the need to comply with stakeholders requests. Most data were made available through body text, but the use of primary databases increased in coincidence with the introduction of complete mitochondrial and next-generation sequencing methods. Our study highlights three important aspects. First, our results imply that researchers' awareness of the importance of openness and transparency for scientific progress may complement stakeholders' policies in achieving very high sharing rates. Second, widespread data sharing does not necessarily coincide with a prevalent use of practices which maximize data findability, accessibility, useability and preservation. A detailed look at the different ways in which data are released can be very useful to detect failures to adopt the best sharing modalities and understand how to correct them. Third and finally, the case of human paleogenetics tells us that a widespread awareness of the importance of Open Science may be important to build reliable scientific practices even in the presence of complex experimental challenges.

Linguistic, geographic and genetic isolation: a collaborative study of Italian populations.

The animal and plant biodiversity of the Italian territory is known to be one of the richest in the Mediterranean basin and Europe as a whole, but does the genetic diversity of extant human populations show a comparable pattern? According to a number of studies, the genetic structure of Italian populations retains the signatures of complex peopling processes which took place from the Paleolithic to modern era. Although the observed patterns highlight a remarkable degree of genetic heterogeneity, they do not, however, take into account an important source of variation. In fact, Italy is home to numerous ethnolinguistic minorities which have yet to be studied systematically. Due to their difference in geographical origin and demographic history, such groups not only signal the cultural and social diversity of our country, but they are also potential contributors to its bio-anthropological heterogeneity. To fill this gap, research groups from four Italian Universities (Bologna, Cagliari, Pisa and Roma Sapienza) started a collaborative study in 2007, which was funded by the Italian Ministry of Education, University and Research and received partial support by the Istituto Italiano di Antropologia. In this paper, we present an account of the results obtained in the course of this initiative. Four case-studies relative to linguistic minorities from the Eastern Alps, Sardinia, Apennines and Southern Italy are first described and discussed, focusing on their micro-evolutionary and anthropological implications. Thereafter, we present the results of a systematic analysis of the relations between linguistic, geographic and genetic isolation. Integrating the data obtained in the course of the long-term study with literature and unpublished results on Italian populations, we show that a combination of linguistic and geographic factors is probably responsible for the presence of the most robust signatures of genetic isolation. Finally, we evaluate the magnitude of the diversity of Italian populations in the European context. The human genetic diversity of our country was found to be greater than observed throughout the continent at short (0-200 km) and intermediate (700-800km) distances, and accounted for most of the highest values of genetic distances observed at all geographic ranges. Interestingly, an important contribution to this pattern comes from the "linguistic islands"( e.g. German speaking groups of Sappada and Luserna from the Eastern Italian Alps), further proof of the importance of considering social and cultural factors when studying human genetic variation.

Demographic histories, isolation and social factors as determinants of the genetic structure of Alpine linguistic groups.

Great European mountain ranges have acted as barriers to gene flow for resident populations since prehistory and have offered a place for the settlement of small, and sometimes culturally diverse, communities. Therefore, the human groups that have settled in these areas are worth exploring as an important potential source of diversity in the genetic structure of European populations. In this study, we present new high resolution data concerning Y chromosomal variation in three distinct Alpine ethno-linguistic groups, Italian, Ladin and German. Combining unpublished and literature data on Y chromosome and mitochondrial variation, we were able to detect different genetic patterns. In fact, within and among population diversity values observed vary across linguistic groups, with German and Italian speakers at the two extremes, and seem to reflect their different demographic histories. Using simulations we inferred that the joint effect of continued genetic isolation and reduced founding group size may explain the apportionment of genetic diversity observed in all groups. Extending the analysis to other continental populations, we observed that the genetic differentiation of Ladins and German speakers from Europeans is comparable or even greater to that observed for well known outliers like Sardinian and Basques. Finally, we found that in south Tyroleans, the social practice of Geschlossener Hof, a hereditary norm which might have favored male dispersal, coincides with a significant intra-group diversity for mtDNA but not for Y chromosome, a genetic pattern which is opposite to those expected among patrilocal populations. Together with previous evidence regarding the possible effects of "local ethnicity" on the genetic structure of German speakers that have settled in the eastern Italian Alps, this finding suggests that taking socio-cultural factors into account together with geographical variables and linguistic diversity may help unveil some yet to be understood aspects of the genetic structure of European populations.

Functional diversity of the glutathione peroxidase gene family among human populations: implications for genetic predisposition to disease and drug response.

AIM: To analyze the human genetic variation of glutathione peroxidases (GPX), estimating the functional differences among human populations and suggesting interethnic differences in predisposition to disease and drug response. MATERIALS & METHODS: Using 1000 Genomes Project data, we analyzed 723 GPX variants in 1092 individuals belonging to 14 populations. Combining functional prediction analyses of coding and noncoding variants, we developed a method to estimate haplotype functionality. RESULTS: GPX rare variants have a higher functional impact than common variants. The frequency among Asian patients of haplotypes associated with normal functionality is significantly higher for GPX1 and lower for GPX3 than for non-Asian patients; no adaptation signals in GPX1 and GPX3 were found in Asian patients. CONCLUSION: GPX1 and GPX3 differences may be associated with alterations in antioxidant capacity and redox regulation, which suggests diverse susceptibility to complex disease and diverse response to relevant drugs in Asians compared with individuals with other ethnic origins.

Detecting genetic isolation in human populations: a study of European language minorities.

The identification of isolation signatures is fundamental to better understand the genetic structure of human populations and to test the relations between cultural factors and genetic variation. However, with current approaches, it is not possible to distinguish between the consequences of long-term isolation and the effects of reduced sample size, selection and differential gene flow. To overcome these limitations, we have integrated the analysis of classical genetic diversity measures with a bayesian method to estimate gene flow and have carried out simulations based on the coalescent. Combining these approaches, we first tested whether the relatively short history of cultural and geographical isolation of four "linguistic islands" of the Eastern Alps (Lessinia, Sauris, Sappada and Timau) had left detectable signatures in their genetic structure. We then compared our findings to previous studies of European population isolates. Finally, we explored the importance of demographic and cultural factors in shaping genetic diversity among the groups under study. A combination of small initial effective size and continued genetic isolation from surrounding populations seems to provide a coherent explanation for the diversity observed among Sauris, Sappada and Timau, which was found to be substantially greater than in other groups of European isolated populations. Simulations of micro-evolutionary scenarios indicate that ethnicity might have been important in increasing genetic diversity among these culturally related and spatially close populations.

Reevaluating a model of gender-biased gene flow among Sub-Saharan Hunter-gatherers and farmers.

In a previous study, we proposed a model for genetic admixture between African hunter-gatherers and food producers, in which we integrated demographic and genetic aspects together with ethnographic knowledge (Destro-Bisol et al. 2004b). In that study it was possible to test the model only using genetic information from widely dispersed and genetically heterogeneous populations. Here we reevaluate the congruence between the model and patterns of genetic variation using an anthropologically and geographically more homogeneous data set that includes Pygmies and farmers from Cameroon, Congo, and the Central African Republic. As implied by the model, the ratios of mtDNA to Y chromosome Nm estimates (effective population size, N, times the migration rate, m; 0.154 in Pygmies and 6.759 in farmers), support an asymmetric gene flow, with a higher Bantu-to-Pygmy gene flow for paternal than for maternal lineages, and vice versa for farmers. Analyses of intra- and interpopulation genetic variation further support the above observation, showing a prevailing effect of genetic drift on maternal lineages and gene flow on paternal lineages among Pygmies, and an opposite pattern among farmers. We also detected differences between patterns for classical and molecular measures of Y chromosome intrapopulation variation, which likely represent signatures of the introgression of Bantu lineages into the gene pool of Pygmy populations. On the whole, our results seem to reflect differences in the demographic history and the degree of patrilocality and polygyny between the two population groups, thus providing further support to our microevolutionary model in an anthropologically coherent framework.

Using forensic microsatellites to decipher the genetic structure of linguistic and geographic isolates: A survey in the eastern Italian Alps.

The study of geographically and/or linguistically isolated populations could represent a potential area of interaction between population and forensic genetics. These investigations may be useful to evaluate the suitability of loci which have been selected using forensic criteria for bio-anthropological studies. At the same time, they give us an opportunity to evaluate the efficiency of forensic tools for parentage testing in groups with peculiar allele frequency profiles. Within the frame of a long-term project concerning Italian linguistic isolates, we studied 15 microsatellite loci (Identifiler kit) comprising the CODIS panel in 11 populations from the north-eastern Italian Alps (Veneto, Trentino and Friuli Venezia Giulia regions). All our analyses of inter-population differentiation highlight the genetic distinctiveness of most Alpine populations comparing them either to each other or with large and non-isolated Italian populations. Interestingly, we brought to light some aspects of population genetic structure which cannot be detected using unilinear polymorphisms. In fact, the analysis of genotypic disequilibrium between loci detected signals of population substructure when all the individuals of Alpine populations are pooled in a single group. Furthermore, despite the relatively low number of loci analyzed, genetic differentiation among Alpine populations was detected at individual level using a Bayesian method to cluster multilocus genotypes. Among the various populations studied, the four linguistic minorities (Fassa Valley, Luserna, Sappada and Sauris) showed the most pronounced diversity and signatures of a peculiar genetic ancestry. Finally, we show that database replacement may affect estimates of probability of paternity even when the local database is replaced by another based on populations which share a common genetic background but which differ in their demographic history. These findings point to the importance of considering the demographic and cultural profile of populations in forensic applications, even in a context of substantial genetic homogeneity such as that of European populations.

Tracing the distribution and evolution of lactase persistence in Southern Europe through the study of the T(-13910) variant.

We investigated the occurrence and intra-allelic variability of the T(-13910) variant located upstream of the lactase gene in 965 individuals from 20 different locations of Italy and Greece. The T(-13910) frequency ranges from 0.072 (Sardinia) to 0.237 (North-East Italy), with a statistically significant difference between North-East Italians and other Italian populations. The comparison of the lactose tolerance predicted by T(-13910) and that assessed by other studies using physiological tests shows a one-way statistically significant discrepancy that could be due to sampling differences. However, the possible role of other genetic factors underlying lactase persistence is worth exploring. The time of the most recent common ancestor and departures from neutrality of the T(-13910) allele were assessed using three microsatellite loci. Time estimates were found to be congruent with the appearance of dairy farming in Southern Europe and the occurrence of a single introgression event. Robust signals of selection can be observed in North-East Italy only. We discuss the possible role of cultural traits and genetic history in determining these observed micro-evolutionary patterns.

Italian isolates today: geographic and linguistic factors shaping human biodiversity.

We briefly review the current status of anthropological and genetic studies of isolated populations and of their micro-evolutionary and biomedical applications, with particular emphasis on European populations. Thereafter, we describe the ongoing collaborative research project "Isolating the Isolates: geographic and cultural factors of human genetic variation" regarding Italian extant geographical and/or linguistic isolates, aimed at overcoming the limitations of previous studies regarding geographical coverage of isolates, number and type of genetic polymorphisms under study and suitability of the experimental design to investigate gene-culture coevolutionary processes. An interdisciplinary sampling approach will make it possible to collect several linguistic isolates and their geographic neighbours from Trentino, Veneto, Friuli, Tuscany, Sardinia and Calabria. This will be coupled with a shared genotyping strategy based on mitochondrial and Y-chromosomal polymorphisms. The results will be analyzed with a focus on the role of geographical and cultural factors in shaping human biodiversity. The aims of the project go beyond the simple reconstruction of the genetic structure and history of the examined groups. In fact, the study will also include an assessment for future bio-medical studies and the development of genetic and bio-demographic databases. Ethical and educational aspects are also foreseen by the project, by using informed consents together with disseminating activities in loco, completed by the creation of a dedicated web site for both scientific and public audiences.

Signal peptide variants that impair secretion of pancreatic secretory trypsin inhibitor (SPINK1) cause autosomal dominant hereditary pancreatitis.

Variants of the SPINK1 gene encoding pancreatic secretory trypsin inhibitor have been described in association with chronic pancreatitis (CP). These alterations are believed to cause a loss of function by either impairing the trypsin inhibitory activity or reducing expression. Here we report two novel SPINK1 variants in exon 1 that affect the secretory signal peptide. The disease-associated c.41T>G (p.L14R) alteration was found in two European families with autosomal dominant hereditary pancreatitis, whereas the c.36G>C (p.L12F) variant was identified as a frequent alteration in subjects of African descent. The functional effects of both alterations and the previously reported c.41T>C (p.L14P) variant were characterized by activity assays and Western blots of wild-type and mutant SPINK1 expressed in human embryonic kidney 293T and Chinese hamster ovary cells. Alterations p.L14R and p.L14P destined the inhibitor for rapid intracellular degradation and thereby abolished SPINK1 secretion, whereas the p.L12F variant showed no detrimental effect. The results provide the first clear experimental demonstration that alterations that markedly reduce SPINK1 expression are associated with classic hereditary pancreatitis. Therefore, these variants should be classified as severe and regarded as disease-causing rather than disease-modifiers.

Phylogeography of the human mitochondrial L1c haplogroup: genetic signatures of the prehistory of Central Africa.

Interindividual variation of human mitochondrial DNA has been extensively studied over the last two decades, and its usefulness for reconstructing evolutionary relationships of extant populations has been proved. However, some mitochondrial lineages still need to be studied using a combination of larger and tailored datasets and increased level of resolution in order to shed light on their origin and on the processes underlying their present distribution. In this study, we analyze the phylogeny of the L1c haplogroup of human mitochondrial DNA using sequence data from hypervariable regions 1 and 2 obtained from 455 individuals (extracted from a total sampling of 2542 individuals) belonging to sub-Saharan African and African-American populations. We propose a substantial revision of L1c phylogeny, by introducing one new sub-haplogroup (L1c4), two new L1c1 clades (L1c1b and L1c1c), and by reassigning the previous L1c1a1 sequences to a clade which we termed L1c5. The new phylogeny encompasses distinct lineages with different evolutionary histories. In fact, based on population frequency, internal variation and mismatch distribution, we propose that L1c1b, L1c1c and L1c2 originated in Bantu ancestors, whereas L1c1a, L1c4 and L1c5 evolved among Western Pygmies. The population structure of L1c is not comparable to any known mitochondrial or, even, Y-chromosomal haplogroup, and challenges the current view that most of mtDNA variation in Pygmies might reflect admixture with Bantu or a persistence of plesiomorphic characters. In fact, the unique feature of the L1c is that it retains a signature of a phase common to the ancestors of the Bantu and Western Pygmies, while encompassing some specific sub-clades which can indicate their divergence. This allowed us to attempt a phylogenetically based assessment of the evolutionary relationships between the two groups. Taking into consideration estimates of the time to the most recent common ancestor of L1c and its clades together with archaeological and paleoclimatological evidence, we propose that the ancestors of Bantu and Western Pygmies separated between 60 and 30 kya.

Brief communication: mtDNA variation in North Cameroon: lack of Asian lineages and implications for back migration from Asia to sub-Saharan Africa.

The hypervariable region-1 and four nucleotide positions (10400, 10873, 12308, and 12705) of the coding region of mitochondrial DNA (mtDNA) were analyzed in 441 individuals belonging to eight populations (Daba, Fali, Fulbe, Mandara, Uldeme, Podokwo, Tali, and Tupuri) from North Cameroon and four populations (Bakaka, Bassa, Bamileke, and Ewondo) from South Cameroon. All mtDNAs were assigned to five haplogroups: three sub-Saharan (L1, L2, and L3), one northern African (U6), and one European (U5). Our results contrast with the observed high frequencies of a Y-chromosome haplogroup of probable Asian origin (R1*-M173) in North Cameroon. As a first step toward a better understanding of the evident discrepancy between mtDNA and Y-chromosome data, we propose two contrasting scenarios. The first one, here termed "migration and asymmetric admixture," implies a back migration from Asia to North Cameroon of a population group carrying the haplotype R1*-M173 at high frequency, and an admixture process restricted to migrant males. The second scenario, on the other hand, temed "divergent drift," implies that modern populations of North Cameroon originated from a small population group which migrated from Asia to Africa and in which, through genetic drift, Y-chromosome haplotype R1*-M173 became predominant, whereas the Asian mtDNA haplogroups were lost.